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Abstract 

The renormalised vacuum expectation values of massless fermion (convention- 
ally call it neutrino) stress-energy tensor are calculated in the static spherically- 
symmetrical wormhole topology. Consider the case when the derivatives of metric 
tensor over radial distance are sufficiently small (in the scale of radius) to justify 
studying a few first orders of quasiclassical (WKB) expansion over derivatives. 
Then we find that violation of the averaged weak energy condition takes place 
irrespectively of the detailed form of metric. This is a necessary condition for the 
neutrino vacuum to be able to support the wormhole geometry. In this respect, 
neutrino vacuum behaves like electromagnetic one and differs from the conformal 
scalar vacuum which does not seem to violate energy conditions for slowly varied 
metric automatically, but requires self-consistent wormhole solution for this. 
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1. Introduction. The possibility of existence of static spherically-symmetrical 
traversible wormhole as topology-nontrivial solution to the Einstein equations has been 
first studied by Morris and Thorne in 1988 They have found that the material 
which threads the wormhole should violate weak energy condition at the throat of the 
wormhole, that is, radial pressure should exceed the density. Moreover, Morris, Thorne 
and Yurtsever |0 have pointed out that averaged weak energy condition (i.e. that inte- 
grated over the radial direction) should also be violated. Since that time much activity 
has been developed in studying the wormhole subject (see, e.g., review by Visser |^) of 
which we consider here the possibility of existence of self- consistent wormhole solutions 
to semiclassical Einstein equations. Checking this possibility requires finding vacuum 
expectation value of the stress-energy tensor as functional of geometry and solving the 
Einstein equations with this tensor as a source. 

Recently self-consistent spherically symmetric wormhole solution has been found nu- 
merically by Hochberg, Popov and Sushkov 0] for the quantised scalar field vacuum 
playing the role of a source for the gravitation. These authors employ vacuum expec- 
tation values of the stress-energy tensor for the scalar field found by Anderson, Hiscock 
and Samuel |^ although Anderson, Hiscock and Taylor have argued (without solving the 
backreaction problem, however) that these values for massive minimally and/or confor- 
mally coupled scalars fail to meet the requirements for maintaining five particular types 
of static spherically symmetric wormholes . As for the experimentally known fields, we 
have calculated in the renormalised stress-energy tensor of electromagnetic vacuum 
with the help of the covariant geodesic point separation method of regularisation [Q. 
It has been found to violate weak energy condition at the wormhole throat in the first 
nonvanishing order in the expansion over the derivatives of metric (the WKB expan- 
sion). This is a necessary condition of existence of self-maintained wormhole solution. 
Important is that this violation takes place irrespectively of the detailed form of metric. 

Also violation of the averaged weak energy condition is necessary (but not sufficient, 
of course) condition for the existence of vacuum self-maintained wormhole. The valid- 
ity of the weak averaged energy condition for the self-consistent solution to the Einstein 
equations has been studied by Flanagan and Wald |P for the massless scalar field. Flana- 
gan and Wald have found that averaged weak energy condition holds for self-consistent 
solution if being additionally averaged transverse to the geodesic using a suitable Plank 
scale smearing function. The Flanagan- Wald result is obtained in the context of pertur- 
bation theory about fiat spacetime, i.e. for the Minkowski topology. In our preceding 
paper |]TD| dealing with electromagnetic field we note that validity of the averaged weak 
energy condition is substantially defined by the space (-time) topology. This condition 
is violated in the first nonvanishing order of WKB expansion over derivatives of metric 
(over radial coordinate) if topology is that of wormhole. 

Important is possibility to have macroscopic wormhole size under proper conditions. 
Flanagan and Wald and Ford and Roman ||1T[ have argued that wormhole size would 
be Plank-scale. These arguments are based on the assumption that coefficient at the 
curvature squared (Weyl term) in the effective gravity action is of Plank scale value. 



Meanwhile, experimental bounds on this coefficient are not very restrictive [|12|, and in 
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pUj] we speculate that this allows the wormhole size to be as large as radius of the Sun. 
A possible large value of the Weyl term might be provided by infra-red contribution into 
effective action from the massless fields, such as electromagnetic one. 

In the given note we perform analogous calculations for the stress-energy of the neu- 
trino vacuum in the wormhole topology. We find that violation of the energy conditions 
takes place in the self-consistent manner analogously to what occurs in the electromag- 
netic vacuum in the wormhole topology 0, |lOl- This is the necessary condition and 
justifies further work towards construction of self-consistent wormhole solution in the 
real physical vacuum of spin 1/2 and 1 fields. 



2. Calculation. The notations we employ are mainly those of ref.|TB[ used in ref.[^. 
The notations for the metric functions r(p), $(p) are read from the following expression 
for the line element: 

ds^ = r\p){de^ + sin^e d(j)'^) + dp'^ - exp (2$)rftl (1) 

The /X, z/, A, ... denote world tensor indices t, p, 6*, while a, /3, 7, ... are purely space such 
indices. We also need notation for the local flat space indices a,b,c,... = 1,2,3. In 
the case of spherical symmetry it is convenient to use decomposition of spinor field 
over spherical waves. Therefore we consider neutrino field as particular vanishing-mass 
case of the four-component Dirac field and take 7-matrices in the local fiat Minkowski 
coordinates x^,x^,x'^,x^ in the usual standard representation: 

(Tq being Pauli matrices, 0"icr2 = io's, ■ ■ ■ ■ We choose tetrad by requiring that 7^ were 
proportional to the expressions which would follow if x^,x^,x^,x^ were treated as usual 
Cartesian coordinates t, x, y, z and t, p, d, (p were treated as usual spherical fiat coordi- 
nates: 



7f = 7oexp ($), 

7p = 7i sin 6 cos + 72 sin 6^ sin -|- 73 cos 0, 

le = (71 cos 9 cos + 72 cos 6' sin — 73 sin 0)r, (3) 
70 = (—71 sin 0-1-72 cos 0)r sin 9. 

The (Tq, are defined by analogous formulas with 7 replaced by a. It proves convenient to 
write bispinor spherical waves in the form 



^ ^ exp(-j./2) ^V+^^ 



Im 



l—m 
21+1 



21+1 ^l,m+l 




exp (— itut). 



l+rn 
21- 



-^rr^i-^^rn+i 



(4) 
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and analogous one with rj^^l^ and ^7+^2+ interchanged. It is imphed that spherical 
function Yim = at m > /. Here / > 1 and the full momentum j = I — 1/2 with m + 1/2 
being it's projection. The Dirac conjugate is ip = ip'^'y'^ and spherical spinors satisfy 

The covariant derivative of a spinor tp-^j, = D^tp = ip^^ — T ^ip can be defined from the 
requirement 

l^,■,u = li,,v - ^\ul\ - ^vl^ + liJ^v = 0. (6) 

This gives 

A = 5i + ^$S7p, ^p = 5p, I^A = 9^ + ^^7A7p, A = Q,cp. (7) 

Description of the formalism is completed by setting the standard expression for the 
action 

S = \\ (V^t'^^^;^ - J~,l^i^)^f^9d^x (8) 
which provides the stress-energy 

% 

T^y = — + V'7i^^;/.) • (9) 



The equation of motion "j^ip-fM = in terms of the two-component function rj = 




takes the form (and the same for ip with ri_Q_ and 1]^Q^ interchanged): 

where we have denoted U = exp(<l>)/r and introduced the new variable z via dz = 
exp {—^)dp. 

Symmetrical separation of the points x, x according to the Christensen's prescription 
1^ for the regularised form of the stress-energy tensor gives, e.g., 



Here {'ip)x, {ip;t)x are the fields ip{x), ip-t{S:), respectively, transported in parallel way from 
X to X along the geodesic. We split the point in radial direction so that x = {t,p,9,(p), 
X = {t,p,e,(p), p- p = 0. 

Now we substitute ip as the field operator expanded in terms of creation and annihila- 
tion operators into (|ll|). As in [§], the operator ordering ip^ . . .ip ^ ■ ■ ■ ■ ■ ■ '^~^) 
is implied. The vacuum expectation values of T^^^ will be denoted as components them- 
selves: this will not lead to any confusion. These values are 
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Anrr 



- exp 



'2" 




(12) 



The tilde denotes the function at p. Here we introduce the three basic sums over quantum 
numbers n (which labels the solutions to the eq.(|TUp) and /.(The summation over m is 
made in closed form thus eliminating the angle 9, dependence and leading to = Tg). 
These sums are 



=Y,i^6i^){v+V+ + V-V-), 

n,l 

Ss = Y.l^e{u){r,^f,.+r^.f,+), (13) 

n,l 

Sa = Y.^^i^)iv+V- - V-V+)- 

n,l 

Here indices n, I at ?7± and 00 are suppressed, 6{uj) = | (l + j^), and rj is normalised as 

-00 

{r]l + r]^_)dz = 1 (14) 



(we use real functions ri±). 

The operator O in ([T0|) is Hermitian w.r.t. the internal product like (|14D and can 
be made self-adjoint by imposing appropriate boundary conditions at z = ±M; then we 
put M — > 00. The two-component eigenvector 7] oi O form the orthogonal set and by 
the completeness property 

Y,r]{z)v'^{z) = l-5{z-z) (15) 

n 

(here 1 is unit 2x2 matrix). Thereby summation over n in the stress-energy is reduced 
to finding the kernel of the operator functions f{0) = 06{0) or 0{O): 

J2f{ujUz)r^+{z) = f{0)6{z-~z). (16) 

n 

These functions can be rewritten as 

9{0) = ^(^l + -^y 0e{0) = \{0 + V0~^). (17) 

Nonlocality upon the action on 5{z — z) arises from the terms with V O"^. Corresponding 
result ai z ^ z can be expressed in terms of Green function G(s, /; z) which satisfies 
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in the following way: 



(18) 



s^ds 



TT 



G{s, /; z, z) 



G{s, /; z, z). 



(19) 



Since the matrix uj'^{z) is diagonal, the WKB expansion for G over the powers of Az = 
z — z and over the derivatives of ^.^(z) does not differ from that obtained for c-number 
LtP'iz) in our previous paper [p!0|]: 



f 1 ^ ii )^\" 
2G{z, z) = exp i-uAz) i - - --^'^ ^ 



2\/2 



UJ 



+Az 



+Az^ 



1 {u^y 1 {ujY ^ 5 



32 uj'^ 

2\/2 



4 Cj3 
4 



;4 

Huj') 



UJ-" 

2\li 



+ 



32 cj6 
32 



-Az"^ 



1 (cu2)" ^ 5 (cu 



2V2 



12 



48 



[00 



2\I2 ' 



32 



(20) 



Here the derivatives are over z. One need to additionally expand over U'^ entering uj'^{z). 

Substitute this into the expression for the stress-energy and introduce instead of s a 
new integration variable q = AzysH^^+lP . Thereby we obtain a collection of terms 
containing powers of derivatives of U over z (that is, of r and $ over p) times coefficients 
of the type 



Az'' 



dq 



f{q) 



UAz qP^q' - U'Az' dqi q^ 



(21) 



where 



/(?) = ?^Z!^6^p(-^0 



2 q ■ 



1=1 



4 sinh^ f 



(22) 



Expansion of the sums (|T^) in powers of Az is achieved by carefully expanding /(g) in 
Taylor series around q = 0. 

When substituting Si, Ss, Sa into ( |T2|) one should also expand f and $ over Az and 
express Az in terms of e = Ap. The result is (the derivatives are over p now): 
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(23) 



(24) 



(25) 



The divergences at e — > can be eliminated by subtracting T^^''^", the stress-energy 
corresponding to the divergent part of the effective action derived by Christensen |Q 
(simultaneously the cosmological constant, Einstein gravity constant and coefficient at 
the Weyl tensor squared C^^xpC^^^^ in the effective action are set equal to their exper- 
imental values). Substituting our metric into the Christensen's formula for the spinor 
field gives: 



3 2 Arrit.c 

in r If' 
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(26) 



(27) 
(28) 



The coefficient at the Weyl term is the only one which is both logarithmically UV and (in 
the considered case of zero mass) IR divergent, and A is the IR cut-off. The Christensen's 
procedure includes also forming half of the sum of the components T^^j corresponding 
to point separations e = ±|e| (above formulas are given for e > 0), thereby only even 
powers of |e| are left (and also |e| rather then e enters the logarithm). Subtracting T^^ 



from T^^s cancels the divergences at e 
we finally obtain 



(this is a useful check of our calculation) , and 
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(30) 
(31) 
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Here L ~ 1; the IR cut off A is fixed only by experiment. The terms of higher orders in 
the derivatives not taken into account in the WKB expansion here can be really omitted 
if the derivatives of $, r are small as compared to unity in the scale of r. 

3. Discussion. The expression for the renormalised neutrino stress-energy tensor 
found differs from the electromagnetic one ||10[ only by the absolute value of the numerical 
coefficients; their signs are the same. Also the difference between the radial pressure 
r = —Tj^'^^'^ and the energy density q = —j'l''^'^'^ at the wormhole throat (where r', $' = 
and r" > 0) is positive. That is, local weak energy condition at the throat is violated. 
Besides that, integrating r — g from p = po < to p = +oo yields the sum of the two 
explicitly positive (in the wormhole topology r'(po) < 0) terms: 



coo 1 j^' 



P=PO 



POO ^'2 



/ (r-f?)exp(-<l>)rfp= -— — -exp(-<l>) 

(32) 

That is, averaged weak energy condition p| is violated as well. The difference from the 
electromagnetic case is in the extent of the weak energy conditions violation which is 
now 16 -7- 20 times weaker. 

The estimate for the wormhole size if logarithm is large compared to unity |10] in the 
presence of A^^i spin 1 and N1/2 spin 1/2 massless fields modifies as 



.2 ^ ^(4iv, + In ( ^ ) . (33) 



Typical size of the wormhole in the neutrino vacuum would be 2 times smaller than that 
in the electromagnetic vacuum. 

It is interesting to compare our expressions for the stress-energy of spin 1/2, 1 massless 
fields with analogous expression for the massless conformal scalar field. As given in refs. 
1^, Q , this expression does not contain the terms of the second order in the derivatives 
of $,r. This seems quite natural because anomalous trace which signals that the stress- 
energy is nonzero in the curved background proves to contain no second order terms 
just in the conformal scalar case. Therefore in our framework we would get t — q = Q 
(the leading terms with no derivatives should vanish m t — q because we must have 
r = at $,r = const due to the t ^ p symmetry in this case) and were to expand to 
the fourth order derivatives whose values are not restricted by the wormhole topology. 
Thus, violation of the weak energy conditions in the wormhole topology with slowly 
varied metric for the scalar vacuum does not take place automatically but should follow 
from the solution of backreaction problem. This is the feature which distinguishes the 
conformal scalar field from the nonzero spin fields. 



References 

[1] Morris Michael S. and Thorne Kip S., Amer.J.Phys. 56 (1988) 395 



8 



Morris Michael S., Thorne Kip S. and Yurt sever Ulvi, Phys. Rev. Lett. 61 (1988) 1446 

Visser M., Lorentzian Wormholes:from Einstein to Hawking (American Institute of 
Physics, Woodbury, 1995) 

Hochberg D., Popov A. and Sushkov S.V., Self- consistent Wormhole Solutions of 
Semiclassical Gravity (Preprint LAEFF 96/25, KSPU-96-03, kr-qc/9701064D to ap- 
pear in Phys. Rev. Lett. 

Anderson P.A., Hiscock W.A. and Samuel D.A., Phys. Rev. D51 (1995) 4337 

Anderson P.A., Hiscock W.A. and Taylor B.E., Stress-energy of a quantized scalar 
field in static wormhole spacetimes (Preprint |gr-qc/9608038| , August 1996), submit- 
ted to Phys. Rev. D 

Khatsymovsky V.M., Phys. Lett. B320 (1994) 234 

Christensen S.M., Phys. Rev. D17 (1978) 946 

Flanagan E.E. and Wald R.M., Phys. Rev. D54 (1996) 6233 

Khatsymovsky V.M., Towards the possibility of a self-maintained vacuum-traversible 
wormhole, to appear in Phys. Lett. B 

Ford L.H. and Roman T.A., Phys. Rev. D53 (1996) 5496 

DeWitt B.S., Phys. Rep. 19C (1975) 295 

DeWitt B.S., The Dynamical Theory of Groups and Fields (Gordon and Breach, 
New York, 1965) 



9 



